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We report the first observation of the charmless two-body mode B —> uK ± decay, and a new 
measurement of the branching fraction for the B ± cjit ± decay. The measured branching fractions 
are B(B ± -> uK = (9.2±l e 3 ± 1.0) x 10~ 6 and B{B ± -> ^tt ± ) = (4.2t^ ± 0.5) x 10~ 6 . We also 
measure the partial rate asymmetry of B — > loK ± decays and obtain Acp = —0.21 ± 0.28 ± 0.03. 
The results are based on a data sample of 29.4 fb _1 collected on the T(4S) resonance by the Belle 
detector at the KEKB e + e~ collider. 
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Charmless hadronic B decays are of interest not only 
for testing our current understanding of heavy quark 
physics, but also as modes to search for direct CP vi- 
olation. The B~ — ► U7r~ and ojK~ decays J|] are dom- 
inated by tree- level and gluonic penguin diagrams 
respectively, illustrated in Fig. [j]. Thus, their branching 
fractions can give us further insight into gluonic penguin 
diagrams, while interference between tree and penguin 
diagrams can lead to a measurable direct CP asymme- 
try. 

In factorization models with N c ~ 2-3, where N c is the 
effective number of colors, the urn mode is larger than the 
uiK mode by a factor of 2 or more Q . This result is borne 
out by further studies in the QCD factorization j|] and 
perturbative QCD (pQCD) § frameworks. In this letter, 
we report measurements of B~ — > ljK~ and W7r~ decays 




u u u ■< II 



FIG. 1: Tree (left) and penguin (right) diagrams for B 
u>K~ and B~ — > lutt~ decays. 



that indicate that the former is more prominent, which 
may suggest the influence of nonfactorized effects. 

The B~ — > ujK~ mode has an interesting history. It 
was first reported by CLEO in 1998 || with 3.9a signif- 
icance, but subsequently superseded by non-observation 
with a larger data set a result that is supported by 
BaBar However, we report here a significant signal 
in this mode. The B~ — > lutt~ mode has been reported 
previously by the CLEO Q and BaBar collabora- 
tions at levels that are somewhat higher than our find- 
ings. The data used in this analysis were collected with 
the Belle detector [|| at KEKB |J , a double storage ring 
that collides 8 GeV electrons and 3.5 GeV positrons with 
a 22 mrad crossing angle. The data sample corresponds 
to an integrated luminosity of 29.4 fb _1 on the T(45') res- 
onance, containing 31.9 million BB pairs, and 2.3 fb _1 
taken 60 MeV below the resonance. 

Belle is a general-purpose detector with a 1.5 T su- 
perconducting solenoid magnet. Charged particle track- 
ing, covering 86% of the total center-of-mass (CM) solid 
angle, is provided by a Silicon Vertex Detector (SVD) 
consisting of three concentric layers of double-sided sili- 
con strip detectors, and a 50-layer Central Drift Cham- 
ber (CDC). Charged hadrons are distinguished by com- 
bining the responses from an array of Silica Aerogel 
Cerenkov Counters (ACC) , a Time of Flight Counter sys- 
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tern (TOF), and dE/dx measurements in the CDC. The 
combined response provides K/n separation of at least 
2.5cr for laboratory momenta up to 3.5 GeV/c. Photons 
and electrons are detected in an array of 8736 CsI(Tl) 
crystals (ECL) located inside the magnetic field and cov- 
ering the entire solid angle of the charged particle track- 
ing system. The 1.5 T magnetic field is returned via a 
flux return that consists of 4.7 cm thick steel plates in- 
terspersed with resistive plate chambers to detect muons 
and Kl mesons (KLM). The Belle detector is described 
in detail elsewhere @]. 

Well reconstructed tracks that are inconsistent with 
being electrons or muons are identified as kaon or pions 
according to a K/ir likelihood ratio (KID), LkI(£-k + 
Ck), where the Ck(-k) are likelihoods derived from the 
responses of the dE/dx, ACC and TOF systems. Candi- 
date 7T° mesons are reconstructed from pairs of photons, 
each consisting of energy clusters greater than 50 MeV in 
the ECL, with m 77 inside a ±3(7 (a = 5.4 MeV/c 2 ) mass 
window around the n° mass jLOj . A mass-constrained 
fit is then performed to improve the 7r° momentum res- 
olution. Candidate ui mesons are formed from 7r + 7r~7r° 
combinations with an invariant mass that is within ±30 
MeV/c 2 of the nominal u> mass |l(J. (The natural width 
of the u> meson is 8.9 MeV.) To futher reduce the large 
combinatorial background from low energy photons and 
7r°s, an co candidate is discarded if the daughter 7r°'s CM 
momenta is below 350 MeV/c. This selection on 7r° CM 
momentum loses 16% of the signal, but removes 60% of 
the combinatorial background. 

We combine an u> candidate with either a K~ or a 7r~ 
track to form a, B~ candidate. As part of this procedure, 
the momenta of the three charged tracks are recalculated 
subject to the constraint that they originate from the 
interaction point. Using the CM beam energy EP^ m = 
\fsj1 = 5.29 GeV and the measured CM energy E™ 
and momentum pg M of the B candidate, we form two 
kinematic variables to select the signal events: the beam- 
constrained mass M bc = \J {E™^) 2 - (p™) 2 an d the 
energy difference AE = E% M - E^ m . 

The major background for this analysis is from contin- 
uum e + e~ — > qq production, where q is a light quark (u, 
d, s, or c). The jet-like continuum events are suppressed 
relative to the more spherical BB events by character- 
ization of the event shape, which is implemented with 
a Fisher discriminant [11 containing six modified Fox- 
Wolfram moments |l2|, 13 1. There are two types of com- 
binatorial backgrounds from continuum events: fake to 
mesons and fake B mesons. The former is suppressed 
using the cross product |P+ xP_ of the momenta of the 
charged pion daughters in the to meson rest frame. The 
latter is suppressed using the B candidate flight direc- 
tion relative to the positron beam axis, and the helicity 
angle of the candidate to meson relative to the B meson. 
The helicity angle, 6>hei, is defined as the angle between 



the B flight direction and the vector perpendicular to the 
lo decay plane in the to rest frame. We use a likelihood 
ratio technique that combines the Fisher discriminant, 
the cross product of the momenta of the charged pions 
from the to, the B flight direction, and the cosine of the 
helicity angle, cos^hci, to suppress the continuum back- 
ground relative to the B — > toh (h = it or K) signal. The 
probability density functions (PDFs) for signal and back- 
ground are constructed using Monte-Carlo (MC) events. 
The background PDFs are in good agreement with those 
determined from on-resonant sideband data (Mb c < 5.26 
GeV/c 2 and \AE\ < 0.3 GeV). With these PDFs, we 
determine signal (£s) and background (£bg) likelihoods 
for each event that are used to form the normalized like- 
lihood ratio 1Z = C§ / (C§ + Cbg)] we discard events with 
1Z < 0.85. This selection retains 50% of the signal while 
rejecting 95% of the continuum background. 

To study background from B decays through the b — ► c 
transition and charmless B decays such as B — > toK* and 
B — ► top, and non-resonant B —> K~ir + ir~ir° decays, we 
used MC samples up to 20 times larger than our data 
sample, assuming the best known branching fraction for 
each decay jl4j. We find negligible backgrounds from 
these decays in the M\> C -AE signal region (Mb c > 5.27 
GeV/c 2 and \AE\ < 0.1 GeV). 

The signal is extracted using Mb c and AE as indepen- 
dent variables in an unbinned maximum likelihood fit for 
events with \AE\ < 0.3 GeV and M hc > 5.2 GeV/c 2 . The 
signal PDF for M^ c is a Gaussian and that for AE is the 
parameterization of Ref . |L5| . The parameters are deter- 
mined from MC simulation and calibrated by the decay 
chain B~ -> D°tt- , D° -> R-ir+ir . The resolutions 
determined from MC are 3 MeV/c 2 for M hc and 24 MeV 
for AE. The PDF of continuum background for M^ c is 
an empirically determined threshold function |l6[ ] that is 
obtained from the sideband data (AE > 0.1 GeV), while 
the PDF for AE is a linear polynomial obtained from the 
data (M bc > 5.27 GeV/c 2 ) before the K cut. PDFs for 
other background sources are included: charmless B de- 
cays that survive the selection criteria and signal events 
with charged kaons misidentified as pions or vice versa. 
(In the last case, the PDFs have the same shape as sig- 
nal except that the central value of AE is shifted by 45 
MeV.) 

The signal yields from the maximum likelihood fit are 
summarized in Table Q. The Mb c and AE distributions 
of candidate events and the best fit curves are shown in 
Fig. §. The signal yields are 18.9±%$ and 10.4^ j events 
for the loK~ and um~ modes, respectively. The expected 
reflection due to ir-K misidentification is 0.7 ± 0.3 (2.0 ± 
0.6) events for the uiK~ (lutt~) mode; the fit gives 0.0±2.4 
(0.0 ± 3.9) events. 

The statistical significance quoted in Table [j] is defined 
as yf— 21n(£(0)/£ max ) where £ max is the maximized like- 
lihood at the nominal signal yield and £(0) is the like- 
lihood with the signal yield fixed at zero. We observe 
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TABLE I: The signal yields, statistical significances (E), effi- 
ciencies (e), branching fractions (B), and the 90% confidence 
level upper limit (UL) of the branching fraction for the ion~ 
mode are listed. The efficiencies include the u) decay branch- 
ing fraction. 
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FIG. 2: The M hc (left) and AE (right) distributions of the 
candidate events (histograms), the best fits (solid curves) and 
signal components (dashed curves). 



18.9 signal events for B~ — > ujK~ with 6.0 a signifi- 
cance and find 10.4 ojt{~ events with 3.3 a significance. 
Since the latter has less than 4 a significance, we use the 
90% confidence level upper limit {Ng L ) of 17.3 events on 
B~ — > ujTt~ yield, determined by integrating the likeli- 
hood as a function of the number of signal events to 90% 
of its total area. 

We study the systematic error associated with the fit 
by varying the parameters in the fitting functions by la 
from their nominal values. The change in the signal yield 
from each variation is added in quadrature to obtain an 
overall systematic error associated with the fit. The sys- 
tematic errors in the detection efficiencies of the uj meson 
and the high-momentum K~ and 7r~ mesons are 8.5% 
and 2.2%, respectively, which are determined from de- 
tailed studies of charged particle tracking, ir° detection, 
and particle identification. A 5% systematic uncertainty 



FIG. 3: The (a) AE distribution and (b) scatter plot of KID 
likelihood ratio versus AE for the uih~ mode. The solid curve 
shows the fit result with the signal components shown by 
dashed curves. 



is assigned to the continuum suppression cut, which is 
obtained by applying a similar procedure to data and 
MC samples of B~ — > D*°ir~ events. The combined un- 
certainty of the efficiency is 10.1%. 

The branching fractions in Table [j] are calculated as- 
suming equal numbers of B + B~ and B°B° pairs in our 
data sample. The uncertainty in the number of BB 
events, 1%, is taken into account and included in the 
systematic error for the branching fraction. The upper 
limit of the branching fraction of ujTt~ decay is calculated 
after increasing Ng L and reducing the efficiency by their 
respective systematic error. 

Our branching fraction result for B~ — > tuK~ is larger 
than that for B~ — ► ujit~ . As a consistency check, we 
also performed the analysis without KID information. 
Figure || shows the AE distribution and a scatter plot 
of the KID likelihood ratio versus AE for the uihT can- 
didates. In these plots, we use the tt~ mass for the high 
momentum hadron track. This causes a —45 MeV differ- 
ence between the peak positions of luK~ and unr~ signals. 
The AE distribution is fitted with loK~ and wtt - signals, 
continuum background, and charmless background com- 
ponents. The signal yields are 17.1 ± 7.7 and 12.1 ± 7.0 
events for ujK ~ and ajn~ , respectively, and are consistent 
with the results using the KID for K/ir separation. The 
scatter plot in Fig. |^(b) shows the distribution of events 
in KID versus AE. The large rectangles, which cover the 
±3cr signal regions in AE and high or low kaon probabil- 
ity, contain enhancements at the appropriate places for 
both modes. 

We also examine the properties of the u> candidates 
to confirm the B~ — > u>K~ signal. The B~ ~ * ojK~ 
signal yield in 7r + 7r~7r° invariant mass bins is shown in 
Fig. [|(a). A clear signal at the uj mass is seen. The 
fitted number of u> mesons is 18.0±5.0 which is consistent 
with the ujK~ signal yield. Figure |^(b) shows the B~ — > 
luK~ signal yield in cos^hei bins. The requirement on 
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FIG. 4: The B~ —> loK~ signal yield in bins of (a) 7r + 7r~7r° 
invariant mass and (b) cosine of the ui helicity angle. The 
solid curve shows the fit result. 



the likelihood ratio has been applied without including 
the helicity angle variable. The distribution is consistent 
with the expected cos 2 $hei distribution. 

We also measure the partial rate asymmetry in B — > 
ujK^ decays to search for direct CP violation. The asym- 
metry is defined as 



A C p = 



N(uK~) - N(ujK+) 
N{ujK-)+N{uoK+)' 



An application of the same event extraction and fitting 
procedure to the B~ and B + candidates separately yields 
7.3±3.5 and 11.2±3.7 events for luK~ and ljK + , respec- 
tively, and an asymmetry value Acp = —0.21 ± 0.28 ± 
0.03. The systematic error includes the uncertainty asso- 
ciated with the fit procedure as well as a contribution of 
1% due to detector bias in reconstruction of positive and 
negative high- momentum kaon tracks. The 90% confi- 
dence level interval —0.70 < Acp < 0.28 is obtained by 
assuming a Gaussian statistical error convolved with the 
systematic error. 

Our combined branching fraction of (13.4^2 9 ± 1-1) X 



10~ 6 for B~ 



jh (h = 7r or K) agrees with CLEO's 



number, (14.3l||±2.0) x 10~ 6 @, although the individ- 
ual branching fractions are not totally consistent. Our 
large B~ — > luK~ branching fraction also disagrees with 
the upper limit of 4 x 10~ 6 reported by the BaBar col- 
laboration Q, although our B~ — > unt~ result is not 
in conflict. We note that BaBar's combined branching 
fraction for B~ — > ujh~ (h — tt or K) is low compared to 
CLEO and our result. 

The large B~ — > luK~ branching fraction and rela- 
tively low B~ — > lutt~ rate cannot be easily accounted for 
either by generalized factorization [2] with iV c ~ 2-3 or 
by calculations based on pQCD j^, I|. To accommodate 
the large B~ — > luK~ branching fraction that we observe, 
it appears that N c has to deviate significantly from 3 jl7) , 
indicating the presence of large non-factorizable effects. 



with the Belle detector, we report the first observa- 
tion of the B~ — » uK~ decay with branching fraction 
B(B- -> luK-) = (9.2t 2 j ± 1.0) x 10~ 6 ; the statistical 
significance of the above signal is 6.0a. We also mea- 
sure B(B- LUTr-) = (4.2+ 2 ;° ± 0.5) x 10~ 6 , with a 
statistical significance of 3.3c The partial rate asym- 



metry for B 



jK ± decays is found to be Acp = 



—0.21 ± 0.28 ± 0.03, corresponding to a 90% confidence 
level interval of -0.70 < Acp < 0.28. 
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